Pauli spin blockade in weakly coupled double quantum dots 
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In a two-level system, constituted by two serially coupled single level quantum dots, coupled 
to external leads we find that the current is suppressed in one direction of biasing caused by a 
fully occupied two-electron triplet state in the interacting region. The efficiency of the current 
suppression is governed by the ratio between the interdot tunnelling rate and the level off-set. In 
the opposite bias direction, the occupation of the two-electron triplet is lifted which allows a larger 
current to flow through the system, where the conductance is provided by transitions between 
one-electron states and two-electron singlet states. Is is also shown that a finite ferromagnetic 
interdot exchange interaction provides an extended range of the current suppression, while an anti- 
ferromagnetic exchange leads to a decreased range of the blockade regime. 
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I. INTRODUCTION 

The non-equilibrium properties of coupled quantum 
dots (QDs) in the Coulomb blockade regime, of both 
serial and parallel configurations^* 2 ^ have been stud- 
ied extensively. Much of the interest in QD systems 
is of fundamental character because of the close anal- 
ogy of QDs as artificial atoms and coupled QDs as 
artificial molecules. 4 Especially, coupled QDs provide 
the opportunity to tune the interdot tunnelling and ca- 
pacitive coupling which change the conductance spec- 
tra and current-voltage (J — V) characteristics of the 
system. QDs displaying ratchet effects have been ex- 
plored both experimentally and theoretically^^ and it 
was recently shown that Coulomb interactions between 
QDs can give rise to ratchet effectsA^ Other phenomena 
have been studied extensively in coupled QD systems 
through the past decade, such as, negative differential 
conductancej2ii2iii I—V asymmetries jiLiS spin-pumping 
and spin-filtering^* quantum phase transitions^* and 
the Kondo effect in serial 15 i 16 i 17 i 18 and parallel 19 i 20 i 21 
coupled QDs. 

Recent transport measurements on double quantum 
dots (DQDs) have revealed technical possibilities to 
utilise the Pauli exclusion principle for the construction of 
efficient current rectifying devices ,22*2* The charge trans- 
port through a series of zero dimensional sites, with bro- 
ken inversion symmetry with respect to the applied bias 
voltage, is suppressed by a combination of the Pauli prin- 
ciple and the Coulomb blockade in one direction while 
permitting it to flow in the opposite direction, thus, 
realising a controllable spin-Coulomb rectifier. In sys- 
tems with absence of inversion symmetry, bias-dependent 
rectification effects can, on the one hand, be provided 
by bias voltage dependent transition rates between the 
states in the system, as was suggested for single molecules 
with appropriately configured molecular orbitals^* and 
for multiple dots suffering from the Coulomb blockade 
effect.' On the other hand, as we show in this paper, 
while the transition rates normally are bias voltage de- 



pendent, this is not a necessary condition to achieve 
the spin-Coulomb, or Pauli spin, blockade and rectify- 
ing properties of the system. Instead, for small inter- 
dot tunnelling rates we find it necessary that the two- 
electron triplet states are aligned (or almost aligned) 
with, at least, the lowest two-electron singlet state, and 
that these states are almost aligned with the lowest one- 
particle states. This condition for the DQD, described 
in more detail later in the paper, agrees very well with 
experimental set-up. While the system displays I — V 
asymmetries and ratchet effects when this condition is 
not satisfied, the origin of those effects is more related to 
normal Coulomb blockade rather than th e more intrigu- 
ing Pauli spin blockade observed in Refs. I'2'2l'23l 

In this paper we undertake a theoretical study of the 
transport properties of a DQD in the regime considered 
in Ref. \22\ . in order to obtain a qualitative picture of 
the origin of the Pauli spin blockade. In excellent agree- 
ment with the experimental report we find that the block- 
ade regime indeed arises because the two-electron triplet 
state in the DQD is occupied with nearly a unit proba- 
bility. The probability for the two-electron triplet state 
to be occupied strongly depends on the ratio between 
the interdot tunnelling rate t and the energy spacing 
Ae = ea — sb, where Ea/b is the single electron level 
in QD^/ B . The DQD system itself provides a range 
of different regimes, of which the one considered here 
is given in the following set up, e.g. |Ae|,J7' » U/2, 
and fi — Taax.{e a, £ b} ~ —U/2, where U' is the interdot 
Coulomb repulsion whereas U is the on-site charging en- 
ergy of QDa/Bi an d n is the equilibrium chemical poten- 
tial. For definiteness, assume that the system is biased 
such that the electron flow is directed according to left 
lead — > QD^ — > QD^ — > right lead. Then, for a small 
ratio 2t/Ae > 0, the probability for an electron in any 
of the two-particle singlet states to exit the DQD to the 
right lead is sufficiently large to prevent a build-up of the 
population in these states. The probability for electrons 
in the left lead to enter the DQD in any of the singlet 
states is of the same order which, however, does not affect 



the inability of charge accumulation in these states. For 
the two-electron triplet states the situation is somewhat 
different, since electrons occupying the triplet has a very 
low probability to exit the DQD to the right lead, whereas 
there is an almost unit probability for an electron in the 
left lead to enter the DQD in the two-electron triplet. 
Hence, there will be an accumulation of electron density 
in the triplet state, which leads to an almost complete 
occupation of this state. The accumulation of electron 
density in the two-electron triplet along with the negli- 
gible probability for occupying the two-electron singlets 
(as well as for the one- and three electrons states) pro- 
vide a blockade of the current, e.g. Pauli spin blockade, 
for a finite range of bias voltages over the DQD. ft is 
important to notice that the occupation probabilities are 
equally distributed among the spin Ms = 1, 0, —1 config- 
urations. Therefore, although one can conclude that the 
Pauli spin blockade is caused by a fully occupied spin 1 
state, it would in general be impossible to read out the 
definite spin state of the DQD in this regime. 

In the opposite case, e.g. small 2tj As < 0, the configu- 
ration of the DQD is mirrored in the sense that the triplet 
state acquires a small probability for occupation since the 
leakage into this state from the left lead is small whereas 
the corresponding leakage to the right lead is large. Thus, 
there cannot be an accumulation of charge in the triplet 
state. In contrast, since the probabilities for electrons to 
enter and exit the DQD via the singlet states are non- 
negligible (as in the former case), the occupation in the 
DQD oscillates between one and two electrons, that is, 
the current is not blockaded by a charge accumulation in 
the DQD. 

The considered DQD system is complicated in the 
sense that we have included on-site charging energy of 
the QDs, the interdot Coulomb and hopping energies, as 
well as the interdot exchange interaction energy. It is 
therefore motivated to describe the physics of the DQD 
in terms of its eigenstate, given in the atomic limit (dis- 
connected to the leads). The current through the DQD is 
thus described by considering single electron transitions 
between the various states. The treatment of the DQD 
dynamics is performed by means of a many-body density 
matrix approach for the population numbers of the eigen- 
states, up to the first order with respect to the couplings 
to the leads. This density matrix approach is similar to 
the one employed by Cota et al~ Also, Coish and Loss 
used this technique in their study of electron spin dynam- 
ics interacting with an environment of nuclear spins via 
hyperfine interactions^ In the model we only take one 
level in each QD, which is sufficient to describe the main 
features of the transport properties for low bias voltages, 
however, in the non-linear transport regime. For larger 
voltages than considered here, it is expected that more 
levels in the QDs participate in the conduction which 
would provide additional understanding to the experi- 
mental observations^ Our model differs from that used 
by Cota et al.f^ in that we also include spin-spin inter- 
actions between electrons in the two QDs. Furthermore, 



while they considered time-dependent transport induced 
by spin-polarised pumping in the DQD, we concentrate 
here on the stationary Pauli spin blockade phenomenon 
which to our knowledge has not been considered before. 

The rest of the paper begins with a definition and dis- 
cussion of the model for the DQD in Sec. [nj an( i the 
employed density matrix approach is described in Sec. 
IIIII In Sec. II VI we analyse the equations for the popula- 
tion numbers, both analytically and numerically in a few 
different cases, while the paper is concluded in Sec. ]V\ 

II. MODEL OF THE SYSTEM 

In order to capture the main features of the DQD stud- 
ied in Ref. I22L we consider two single-level quantum dots 
(QDs) with the intradot Coulomb energies U A / B , inter- 
acting via the hopping interaction t, and interdot charg- 
ing energy U' a nd exchan ge J. Thus, we introduce the 
model, c.f. Refs. Il3l27l28l 

"Hdqd = ^ £Acrd Aa d Acr + U A n A ^n Al 

a 

+ ^ £ Bcyd} B(y d B<y + U B n B ]n B i 

a 

+ (W - J/2) (rut + n Ai ){n B] + n Bl ) 
-2Js A -s B + Y,t(d Aa d Ba + H.c), (1) 

a 

where d Aa , Bu (d Aa / Ba ) creates (annihilates) an elec- 
tron in QD A / B at the energy e Aa / Ba and spin pro- 
jection a =T>4) whereas n Aa — d Aa d Aa and likewise 
for n Ba . The sixth term in Eq. Q accounts for 
the spin-spin interactions between the QDs with = 
(1/2) J2 aa > d Aa a aa 'dAa', and analogous for s B , where 
a = (a x ,a v ,a z ) is the vector of the Pauli spin matri- 
ces. Below we will only consider spin-degenerate QDs, 
hence, we set e Aa = e A and e Ba = e B . This assumption 
is consistent with the experimental set-up of Ref. [2^ since 
there were no magnetic field applied or (any significant) 
spin-orbit splitting of the QD levels. 

Because of the inherent complexity of this model it is 
motivated to transform the Hamiltonian TLdqd into di- 
agonal form. Therefore, we need to find the eigenstates of 
the system, Eq. which are given on the form \N, n), 
representing the nth state of the TV-electron configura- 
tion, see Table [I] Obviously, there are 16 different states 
of the DQD. The corresponding energies are denoted by 
E 01 = 0, Ei ni E 2n , E 3n , En, respectively. Thus, the 
model given in Eq. Q can be written in diagonal form 
according to 

4 

Ti-DQD = ^J-SiVn^AT- (2) 

N=0 n 

Here, we have introduced the Hubbard operators22i3i 

X-NN' ~ l-^i n ) n '\ an d the diagonal transitions h 7 ^ = 



2 



X™ N . The diagonalisation process becomes clear by con- 
sidering the expansion 

d A * = ^^|iV,n)(iV ) n|^ CT |iV> / )<iV , ,n'| = 

NN' nn' 

— E y^(^a-)jVJV+1^7VW+li 
N nn' 

where we have put (cUoOjvlv+i = (N,n\d a \N + 1, n'), 
for a shorter notation. A similar transformation can be 
made for the operator dsa- The summation runs over 
all states \N,n) and \N + l,n') in the DQD, where we 

have used that (dA*)™^' = ^JV'iv+iC^^^+i- Tne ex " 
pansion above along with the algebra of the Hubbard 

operators, i.e. X^ n N ,X^ M , = 5 N >M8n> P X N P M' y ield tne 
desired transformation of the Hamiltonian to the form in 
Eq. (J2J). It should be noticed that the transition matrix 
elements (N, n\d,Aa/Ba\N' , n') are crucial for the trans- 
port properties of the system, as will become clear later 
in the paper. 

The total system, constituted of the DQD and the con- 
tacts to the left (L) and right (i?) of the DQD is given 
by 

H = E e ka c\ a c k(J + Hdqd 

+ E ( V kv c lcrdAa + H.C.) 
ka£L 

+ E ( v kvcLdBa + H.c.) (3) 
kaGR 

where c\ a (cfc ff ), ka € L/R, creates (annihilates) an elec- 
tron in the left /right lead at the energy Ska- The last 



TABLE I: Schematic representation of the eigenstates for 
the model Hamiltonian TLdqd in Eq. Q. Here, \a) = 
4-|0, 1), where jO, 1} = |0> |0> is the empty state of the DQD, 
whereas \a)\a') = \a) A \a') B = d^d^JO, 1), and | Tl)|0) = 
d^id^|0, 1} and correspondingly for QDs. The two-particle 

states |*£f ) = [| T)| I) - I 1)1 WV2, \^ =0 ) = | U>|0>, and 
|$f=o) = 1 0> | U). The coefficients a„, p„, A n , B n , C„, «„, 
and X„ depend on the internal parameters e Aa jBa, U, U' , J 
and t of the DQD. 



N= 1 


|l,n)=a n |t>|0>+/3»|0>|t), n = l,3 
|l,n}=a„|l}|0)+/3 n |0}|t), n = 2,4 


N — 2 


|2,1> = |T)|T>, |2,2> = |4>|4) 

|2,3> = [| T>| 1> + Ii)l T)]V2 

\2,n) = A n \^ ) + B n \^ =0 ) + C„|$f =0 ) 

n = 4,5,6 


N = 3 


3,n) = Kn\ tl) T) + A„| t) Tl), n = l,3 
|3,n> =Kn\ Tl) 4) + A„| I)| Tl), n = 2,4 


N = 4 


IT1) = 1 Tl) Tl) 



two terms express the tunnelling between the contacts 
and the DQD. In order to obtain a uniform description 
of the physics in the system, it is desired to express these 
in terms of Hubbard operators as well. By the above ob- 
servations, we note that the tunnelling can be written as 

H T = 5>fcrcL E J2( d ^ n N + l X NN + l + H.C.] (4) 
ka N nn' 

where d a = d A a or d„ — dsa depending on whether 
ka G L or ka G R. Note that the spin indices on Cka and 
d a have to coincide, since we neglect possible spin-flip 
transitions in the tunnel barriers. Such transitions, how- 
ever, may occur within the DQD but will be neglected 
here. 



III. DENSITY MATRIX FOR THE 
MANY-BODY OCCUPATION NUMBERS 

In the model we consider, there are 16 different many- 
body states, cf. Table [I] and hence there are 16 occupa- 
tion numbers (P/v™, N — 0, ... ,4) to calculate. These 
occupation numbers provide the probability for the state 
\N, n) to be populated. Here, this will be done by means 
of a density matrix for the many-body occupation num- 
bers, e.g. p = {X™ N ,}Nn.N'n'- However, since we are 
only interested in the evolution of the occupation num- 
bers up to second order in the hybridisation v^a , we can 
neglect the influence from the off-diagonal components 
in the density matrix. We can safely neglect all the off- 
diagonal terms, in the second order approximation, since 
i) the off-diagonal terms only give contributions to the 
diagonal components of the fourth order and higher in 
Vka, and ii) the off-diagonal components are multiplied 
by off-diagonal transition matrix elements (to the fourth 
order) which in general are very small compared to the 
diagonal transition matrix elements. 

In the present case we only consider stationary cur- 
rents through the system, hence it is sufficient to in- 
vestigate dp/dt = subject to the boundary condition 
1 = J^Nn-^Nn- The equations take the form 

ka n' 

— E WfeCT (^ CT )AW+l( C fea-"^V7V+l) J i 
n' ' 

J^4i = -^ImEE^^^aKcL^" 1 )- (5) 

ka n 

In the second equation above the electron number N = 
1, 2, 3. To the first order of approximation (in the Marko- 
vian limit), which is sufficient for describing the current 
in the second order approximation, the imaginary part 



3 



of the correlation functions (ct-Xj^,), N' — N ± 1, are 
given by 

Im(cL^) = -™tM)iUfi{£k«)PNn (6) 

where f^{e k a) = fi^ka - fi x ), k <7 e x = L,R, 
is the Fermi function whereas /"(fffeo-) = 1 — 



9 P 
dt PNn 



IV. ANALYSIS OF THE POPULATION 
NUMBERS 

The main task of this paper is to show that the Pauli 
spin blockade is caused by a full occupation of the two- 
electron triplet in the DQD. For this purpose, we tune 
the system into the conditions that were established in 
the experimental situation of Ref. l22l with two identical 
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FIG. 1: Equilibrium distribution of the one-electron tran- 
sition energies A Nn N / n i in terms of Ae, for ea = —U/2, 
Ae = ea -e B = U/2, U' = U/2, 2t/Ae < 1, and J = 0. 
In the absence of magnetic field we have Ai2,oi = Au.oi, 

Al4,01 = Al3,01, A26,12 = A26.ll, and A26,14 = A26,13- 



fx( e ka), and A N > n ^ Nn = E N > n > - E Nn denote the en- 
ergy for the transition XJjnyy. Introducing the cou- 
plings to the left and right leads by ^^^N'ri' = 

2 71 " ^2kaeL/R \ v ko (^(t)jvjV' | 2< ^(AjV'n',iV>i ~ £fecr), w< 3 hnd 

that the equations for the many-body population num- 
bers can be written as 



(7a) 



(7b) 
(7c) 



QDs. There, the internal parameters of the DQD were 
configured such that the levels of the two QDs were sepa- 
rated by Ae = ea — £b = U/2, where U — Ua — Ub, and 
that £a~ M = —U/2. Moreover, the interdot Coulomb re- 
pulsion U' = U/2 and the interdot hopping 2t/Ae <C 1. 
We note that the exchange interaction parameter J in 
the experiment was estimated to be much less than any 
other parameter in the system, e.g. J <C ksT ~ U/10, 
thus we set J = 0. The one-electron transition energies 
between the different states are graphically displayed in 
Fig. ^ Under those circumstances we find that only 
|A 2n> i n > ~ Ml U/2 for n = 1, . . . , 5, and n' — 1, 2, and 
|^3n,26 _ A*l ^ U/2 for n — 3,4, where \i is the equi- 
librium chemical potential of the system. For all other 
transitions we have |Ajv nj jv_i n / — fi\ > U/2, c.f. Fig^ 



A. The Pauli spin blockade regime 

Now, we let the bias voltage^ 3 . eV = (1l — M-R € 
[0.1, 1]U such that /+(A 2 „,!„0 = 1, /^(A 2n , ln = 

0, n = 1,...,5, n' = 1,2. In this regime, we also 
note that /+(A ln ',oi) = 1 for all n', /+(A 2 „,i n /) = 

1, n = 1,...,5, n 1 = 3,4, /+(A 26 ,in') = for all 
n', /+(A 3n /, 2n ) = 0, n' = 1,...,4, n = 1,...,5, 
/+(A 3 „,, 26 ) = 1, ri = 1,2, #(A 3n ,, 26 ) = f R (A 3n , i26 ) = 
1, n' = 3,4, and that /+(A 4 i, 3n ') = for all n'. Then, 
from Eq. (|7a|l and l|7c|l it follows that P01 = P41 = 0, as 
expected since the transition energies Ai„/.oi and A 4 i )3 „' 
lie far from resonance. It also follows that P 2 6 van- 
ishes, which is clear from the fact that the transition 



J 



— J! r 01,lnl/x (Aln.Ol)-Poi ~ f x (Al„,0l)-Pln] 

Xn 

= ^^\^ T N-ln',Nnifx^ Nn ' N ^ ln '^ PN - ln ' ~ fx ( &Nn,N- In' )PjVn] 
X n' 

~ ^ rjVn,JV+ln' [/>^(AjV+ln',A'n)-PA r n ~ f% (AjV+ln',iVn)-PjV+ln'] ) 
n' ' 

= ^E r 3n,4l[/ x + (A 4 l,3„)P3n - /" (A 41 , 3n )P 41 ] • 

X™ 

I 
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energies A 2 6,in' — /i > U/2. Therefore, there cannot oc- 
cur any transitions between the state |2, 6) and any of 
the three-electron states, although the transition ener- 
gies hr < A 3n / : 26 < Hl, n = 3,4. Taking this obser- 
vation together with the fact that A^n/^n — H > U/2, 
n = 1, . . . , 5, n' = 1, . . . , 4, it is clear that the population 
numbers P^ n i — for all n' . This implies that the pop- 
ulation numbers Pi n > of the one-electron states \l,n'), 
n' = 3,4 vanish. Hence, the only non-vanishing popula- 
tion numbers are P\ n ', n 1 = 1,2 and P^n, n = 1, . . . , 5, 
which lead to the equations 

D Z^n=l 1 ln',2n r 2n / 1 o (Q ^ 

An' = -=b ^ . n = 1 -2 J (N 

In' ,2m 

n'=l,2 

n=l,...,5, (8b) 

where we have used that f//(^2n.in') = 1 an d 
/fl(A2 ra ,in') = 0. The three configurations of the two- 
electron triplet, e.g. |2, ti), n = 1,2,3, have to have 
equal probability (spin degenerate system), hence, we set 
Pin = Nt/3, n = 1, 2, 3. The transition matrix elements 
(d a )i 2 n are in general given in terms of the internal pa- 
rameters of the DQD, see Table [H] Using these we find 
that P\x — Pi 2 = Ni/2, which is expected since the DQD 
is spin degenerate. From Eq. JHJ we then also find that 

2 T R ( a \ 2 

N, - 5fr y N T (9.) 

p - = Kfe'i) 2wT ' n=4 ' 5 - (9b) 

where L„ = /3j4„/ %/2 + aB n and P„ = aA n /\[2 + (3C n , 
which are bounded and finite for all £ = 2t/ Ae and slowly 
varying functions of £. Thus, by the normalisation condi- 
tion, now reduced to 1 = JVi + Nt + P24 + P25, we finally 
arrive at 




TABLE II: Matrix elements for the transitions X™ 2 n given in 
terms of the parameters of the eigenstates of the model, Eq. 
0, c.f. TableQ Here a = a n , = f3 n , n = 1, 2. 
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FIG. 2: (Colour online) The two-electron population numbers 
N T (solid), P24 (dotted), and P25 (dashed), for £ = 0.1 (bold) 
and £ = 0.5 (faint), holding Ae constant, and k B T ~ [7/10. 
Other parameters as in Fig. Q 

Hence, the probability for populating the two-electron 
triplet in the restricted bias voltage interval can be given 
in terms of the internal parameters of the DQD and the 
couplings to the leads. Especially, we observe that Nt 
approaches 1 for £ < 1, since a 2 = £ 2 /[(l + y/l + £ 2 ) 2 + 

e\ « Oand /3 2 = (l + v / T+F) 2 /[(l + V / l+^) 2 +^ 2 ) « 
1. This is consistent with Eq. |J3J, which shows that 
iVi and Pjn, n = 4,5, tend to zero in the same limit. 
The main conclusion then is that the occupation of the 
two-electron triplet strongly depends on the ratio be- 
tween the interdot hopping strength t and the relative 
level separation Ae, e.g. £, and that the probability of 
a fully occupied two-electron triplet approaches unity as 
£ — > 0. One should note, however, that this conclusion 
follows independently of a possible bias dependence of 
the transition matrix elements. Including such a depen- 
dence into the consideration only strengthens our conclu- 
sion. A similar analysis can be performed with a finite 
exchange interaction parameter J, which will result in 
an equivalent conclusion of the existence of a Pauli spin 
blockade regime, see Sec. IIV Dl In Fig. |3we display the 
bias dependence of the two-electron population numbers 
Nt, P2n, n — 4,5, for two different values of £, and it 
is clear that the two-electron triplet state is almost fully 
populated in the considered regime of bias voltages for 
small £, whereas it becomes depopulated for increasing 
£. It is also clear from the figure that P2„, n = 4, 5, are 
negligible in the spin-Coulomb regime. 

The current I through the DQD will be suppressed by 
the high population of the two-electron triplet state. In 
the stationary regime I = (It, — Ir)/2, II = —Ir, where 
the current I x flowing from the x lead to the DQD is 
given by 

h = -|tr Im j rx[4+HG>H + f-(u)G<(u)]<kj, 

(11) 

where the greater and lesser Green functions of the 
DQD in the zeroth order approximation are given 
b y 10 ' n i 32 G^ n N+ln ,(tu) = i2TrP N+ln 'S(oj - A N+ln ^ Nn ) 
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bias voltage (eV/U) 



FIG. 3: (Colour online) The I — V characteristics of the DQD 
system for various values of £, holding Ae constant, fesT ~ 
U/10, and r = 2F L . The lower right inset shows the current in 
the Pauli spin blockade regime. The upper left inset shows the 
current for bias voltage dependent transition matrix elements 
imposing a 10 % voltage drop between the QDs, and £ = 0.05 
(same scale as main panel). 

and G^ n N+ln ,(u) = -i2irP Nn 5(uj-A N+ i n ^ Nn ), respec- 
tively. Here, the coupling matrix T a includes the transi- 
tion matrix elements of the DQD. Using the above anal- 
ysis for the population numbers we find that the current 
Ij^ in the Pauli spin blockade regime reduces to 

I? = |r»(j-)V/2 + L; + ^]iV T (12) 

where the last line is obtained in the case of symmet- 
ric couplings, T L = T R = T/2. As expected, the cur- 
rent is proportional to the population number Nt of the 
two-electron triplet state. More importantly, though, is 
that the current strongly depends on the ratio (a//?) 2 = 
£ 2 /(l + yl + £ 2 ) 2 , showing a strong suppression of the 
current for weakly coupled QDs (£ <C 1). 

In Fig. El we show the I — V characteristics of the DQD 
system in the given configuration for different £, holding 
Ae constant. This shows a clearly suppressed current 
in the forward biased system, e.g. eV = (Al — A*r > 0. 
The lower right inset shows the forward current in the 
Pauli spin blockade regime on a logarithmic scale, and 
for small £ it is readily seen that the current is about 1 
— 4 orders of magnitude smaller in the blockaded regime 
than in the linear regime around zero voltage. This ver- 
ifies the experimental observation presented in Ref. 12^ . 
since bias voltage dependent transition matrix elements 
merely provide a slight modification of the current, see 
upper left inset of Fig. El where we have followed the 
experimental estimate of a 10 % voltage drop between 
the QDs. The bias dependence of the transition matrix 
elements is provided by assuming a bias dependent shift 
of the single particle levels in the QDs. 

The current I~£ is non-vanishing in the spin-Coulomb 
regime because of the finite hopping interaction t be- 
tween the QDs, which causes the triplet state population 
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FIG. 4: a) Differential conductance (dl/dV) of the DQD for 
varying positions of the chemical potential /i£ [— 1, 2] U. The 
areas between the grey (black) ridges correspond to dl /dV = 
0. b) Contour plot of the population number Nt = 1 as 
function of the bias voltage and the position of the chemical 
potential. The truncated diamond in the bias voltage range 
[0, 1.5] U marks the region with Nt ~ 1. Here, £ = 0.05 while 
other parameters as in Fig. |5] 



number Nt to be slightly less than unity. The remain- 
der of the available population is distributed among the 
one-particle states |l,n), n = 1,2, and the two-electron 
singlets |2,n), n = 4, 5. The non-vanishing population 
number N± implies that there is a finite probability for 
one electron to exit the DQD to the right lead, and there- 
after another electron enters the DQD from the left. The 
leakage current in the blockaded regime grows quadrati- 
cally with £ since (a//3) 2 = £7(1 + y/l + £ 2 ) 2 , while Nt 
decays by 1/(1 + (a//3) 2 /3), see Figs. and El 

In the above analysis of the population numbers, the 
equilibrium chemical potential \x = 0. Shifting the po- 
sition of fi corresponds to applying a gate voltage over 
the system, which changes the positions of the QD lev- 
els relatively /i. In doing this we find that the Pauli 
spin blockade regime extends over a range of gate volt- 
ages, where Nt ~ 1 and the current through the DQD 
is small. In Fig. El b) we show a contour plot of Nt 
as function of the bias voltage and n, readily showing 
the expected domain of the Pauli spin blockade regime 
(truncated diamond in the bias voltage range (0, 1.5)U). 
The corresponding differential conductance of the DQD 
system is plotted in Fig. El a), where the grey and white 
ridges correspond to dl/dV > while the black ridges 
corresponds to dl/dV < 0. The areas between the ridges 
correspond to vanishing differential conductance. 

The population number Nt of the two-electron triplet 
state, which strongly depends on (a//3) 2 , depends on the 
couplings T L/R to the leads as r fl /r L , see Eq. iJTUJ. 
Hence, assuming proportional couplings, i.e. r^/r L = 
7 > 0, such that 7 ~ (/3/a) 2 , provides a lifting of the 
Pauli spin blockade even for weakly coupled QDs. Ac- 
cordingly, for £ — > one finds that Nt — * 3/5. In Fig. 
0we plot the population number Nt (bold) for different 
asymmetries of the couplings to the left and right leads, 
and it is readily seen that Nt decreases as 7 is increased. 
Especially, for 7 = (j3/a) 2 it is seen that Nt is close to 
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Nx, Nt, P25 in terms of P24, which give 
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FIG. 5: (Colour online) The two-electron population number 
Nt (bold) for different asymmetries 7 = p fl yp L Q f the cou- 
plings to the leads. The sum of the one-electron population 
numbers Nx (faint) are plotted for reference. Here, £ = 0.05, 
whereas other parameters as in Fig. [5] 



3/5. The population numbers p2 n , n = 4, 5, remain neg- 
ligible as 7 is increased since they do not depend on this 
ratio. The remaining population (1 — Nt) is distributed 
among the one-electron states \ln), n — 1,2, see Eq. 
(|9al) . since Nx oc jNt- In Fig. 0we plot the population 
number Nx (faint) of the one-electron states, illustrating 
that the lifted Pauli spin blockade results in an increasing 
population in the one-electron states. 

The absolute magnitude of the current, however, is 
not increased by the lifting of the Pauli spin blockade. 
In the limit 7 = x(f3/a) 2 and ( « 1, where x > is 
an arbitrary constant, the current becomes (using T — 

r L + r R = (1 + 7)r L ) 



P 



4xe 



8ttx 



1 eT 



(3 + 2x)h 



3 + 2x 1 + 7 h ' 



(13) 



showing that the current remains small (decreasing) 
for increasing asymmetries 7 of the couplings to the 
leads. Nevertheless, as we will see in the next para- 
graph, Sec 



IV Bl the current /£" in the reverse bias 
direction (in the range — [1,0.1] U) becomes in this limit 
w [Tr/(l+j)]-eT/h. It is, thus, clear that the ratio be- 
tween the forward and backward currents, e.g. P^V-^zT 
is a bounded and monotonically increasing function of 
the parameter 7. Hence, although the absolute magni- 
tude of the current is not increased by the lifting of the 
Pauli spin blockade caused by the asymmetric couplings 
to the leads, the forward current increases relatively the 
backward current for growing asymmetries T R /T L . 



B. Reverse bias 

An analogous investigation of the case eV = hl — Hr G 
— [1,0.1]{7 shows that the population number Nt of the 
two-electron triplet approaches zero for ( < 1. This 
is easiest seen by expressing the population numbers 



(14a) 
(14b) 
(14c) 



These equations together with the normalisation yields 



P 2 4 = U + T=i 



P14 



T L 3 fa 



R5 

L 5 



(15) 

Thus, for small £ it is readily seen that Nt ~ 0. This also 
shows that the population of the one- and two-electron 
states I l,n'), n' = 1,2 and |2, n), n = 4, 5, respectively, 
are finite, meaning that the electron population in the 
DQD varies between 1 and 2. The reverse bias case 
corresponds to negative biases in Fig. showing the 
small triplet state population number Nt and finite sin- 
glet population numbers P24 (dotted) and P25 (dashed). 
In the considered regime, the sum P24 + P25 < 1, and 
since Nt is negligible the remaining population is dis- 
tributed among the one-electron states, e.g. population 
number Ni. Hence, the electron population in the DQD 
flips between 1 and 2 which allows a finite current to 
run through the system. The current T£ in the reverse 
direction, given by 



2e 



T L Ll[l + (R 5 /Ri) 2 }P, 



24 ■ 



eF 



2ir^-Ll[l + (R 5 /R 4 ) 2 ]P 2 4. 



(16) 
(17) 



is expected to be substantially larger than in the forward 
direction. Here, the last line is again obtained in the case 
of symmetric coupling. 



C. Linear regime 

In the linear regime, e.g. eV £ (— 0.1, Q.1)U, we can 
employ the equilibrium result 

N ± = 2e (A24 - 11 ^ l)/fcBT P 24 (18a) 

N T = 3 e ( AM .n-A 21 , nVfceTp^ ( 18b ) 

P25 = e (A24 ' 11 - A25 ' ll)/fesT P 2 4, (18c) 

which gives 

P24 = {1 + 2e*- A24 ' 11_Al ^ fcBT + 3 e ( A24 ' 11 ~ A21 - 11 )/ fc - EiT 

_|_ e (A 2 4,ii-A25,u)/fcsTj-l^ rtgj 

Since A 2 4 y xi < M = 0, see Fig. ^ we see from these 
relations that P24 ~ 1 whenever |A 2 4,n| 3> kgT which 
is satisfied for £ > 0.1 and ksT ~ U/10. Accordingly, 
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the population numbers in Eq. i|18|) wm be exponen- 
tially suppressed by the same condition on £, since the 
gap between the energy of the lowest two-electron singlet, 
e.g. 1 2, 4), and the other two-electron energies increase 
roughly by £ 2 . Hence, for £ > 0.1 the current will be 
blockaded by the lowest two-electron singlet in a finite 
bias voltage interval around equilibrium. The plots in 
Fig. |21 and |3] illustrates the blockade because of the fully 
occupied singlet state P24 in the region around zero bias 
voltage for large £. It is also clear that P 2 4 — ► 1/7 as 
£ — ► 0, since the energy distance between A 2 4,n and 
fi, A2n,ii! ti — 1,2,3,5, approach zero. In the limit of 
weakly coupled QDs, e.g. £ — > 0, the equilibrium popula- 
tion numbers satisfy Ni = 2/7, Nt = 3/7, P24 = P25 = 
1/7. The conductance of the system in the linear regime 
for weakly coupled QDs is therefore substantially larger 
than in the Pauli spin blockade regime. 



D. The role of the interdot exchange 

In the previous sections we have considered the ex- 
change interaction parameter |J| <C kgT, in which case 
the singlet-triplet splitting is negligible. Here, we will 
consider the case |J| > 2ksF ', where J > gives a 
ferromagnetic interdot exchange in our model whereas 
J < provides an anti-ferromagnetic interdot exchange. 
The energy of the triplet state is given by i? 2 „ — 
Ea + £b + U' — J, n = 1,2,3, whereas the lowest sin- 
glet state energy is roughly E s = min n= 4 ) 5, 6 {-E l s'2n} ~ 
£a + £b + U' + J — F(t, Ae), where F(t, As) is a positive 
function. For weakly coupled QDs, e.g. £ <C 1, we find 
that F(t, Ae) w £. Hence, in the configuration of the 
DQD system of interest in this paper and for low tem- 
peratures, the lowest singlet state is the ground state of 
the DQD whenever J < £/2 while the triplet becomes 
the ground state for larger J. 

Here, the task is to demonstrate, for £ < 2/cgT, that 
there is a finite range of bias voltages population prob- 
ability of the triplet state is almost unity, at least for 
interdot exchange interactions —Q.5(U — U')/2 < J < 
(U — U')/2. Indeed, when the system is biased such that 
fi R < A 2 „,i n / < hl, n — 1, . . . , 5, n' = 1,2, while all 
other transition energies lie above or below both fj.^ and 
fiR, we can apply the same arguments as in Sec. IIV Al 
Then, it follows that Nt — > 1 whenever £ — > 0. This 
argument is valid for all 2J/(U - U') S (-0.5, 1). It is 
clear, however, that this range of interdot exchange in- 
teractions contains more interesting features. 

First, consider the range < J < £. While the lowest 
singlet state remains the ground state of the DQD in 
this regime, the transport properties of the system are 
not significantly changed from the case of vanishing J. 
Hence, the Pauli spin blockade regime is again found for 
positive biases, however, slightly extended with larger 
J since the transition energies A2„,i n ' — fi < 0, n = 
1,2,3, ri = 1,2. 

Secondly, for J well above £, the triplet state acquires 
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FIG. 6: (Colour online) I — V characteristics a) and triplet 
state population probability b) for various strengths of the 
interdot exchange interaction. Here £ = 0.05 and other pa- 
rameters as in Fig. [5] 



the lowest energy among the two-electron states, e.g. 
A2„,i n ' — m i n m=4,5,6{A 2m ,i„'}j n — 1,2,3, n = 1,2. 
Hence, depending on the value of J in the interval 
(£, [U—U']/2), the Pauli spin blockade regime is extended 
to include also equilibrium as well as a finite interval of 
negative biases. The reason for this is to be found in 
the matrix elements of the transitions between the one- 
electron states and the lowest singlet states, while the 
matrix elements for transitions between the one-electron 
states and the triplet are independent of J. For small J, 
the two lowest singlet states are almost equally weighted 
on the Fock states \^s=o) an d l*f=o) ( see Table[j]for the 
definitions of these states). Hence, those state are cou- 
pled to the left and right lead with almost equal strength. 
For large J, however, the weights of the singlets are re- 
distributed among the Fock states such that the lowest 
singlet acquires a larger weight on QD^, e.g. on the state 
|$g =0 ), whereas the other is more strongly weighted on 
l^s^o) ■ This, in turn, leads to that electrons tunnelling 
to the lowest singlet state in the DQD from the left/right 
lead is small/large which prevents accumulation of elec- 
tron density in this state for biases such that A 2ni i„' — 
fi x < 0, n = 1,2,3, n' — 1,2, and min x= L,fl{/i x } < 
min m=4 , 5i6 {A 2mi i„/} < max x=Lji? {/i x }. For larger neg- 
ative biases, e.g. such that /x/, < A2 n ,in' < M-Ri n 
1.2..'i. n! = 1,2, the population in the triplet is lifted 
and the current can flow through the DQD via transitions 
between the one-electron states and the lowest singlet. 

Third, a anti-ferromagnetic exchange, e.g. J < 0, 
yields a smaller range of positive bias voltages (com- 
pared to the case J = 0) where the triplet state is 
fully populated. This is understood as an effect of that 
A 2nj i„' — fi > 0, n = 1, 2, 3, n' = 1,2. Hence, a larger 
bias voltage is required (than in the case of J = 0) to 
enter the Pauli spin blockade regime. 

In Fig. we display the I — V characteristics a) and 
triplet state population probability for several values of 
the interdot exchange interaction parameter J. The plots 
in the figure clearly shows that the transport proper- 
ties of the system for \2J/(U — U')\ < 0.1 resembles 
the case of negligible J. It is however, also clear from 
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these plots that the Pauli spin blockade regime is en- 
tered for lower biases as J > and for higher biases as 
J < 0, as expected from the above discussion. The case 
2J/(U — U') = —1 has been included for completeness, 
where the blockade around equilibrium is caused by a 
full occupation in the lowest singlet state, that is, con- 
ventional Coulomb blockade. 



V. CONCLUSIONS 

We have theoretically verified the experimentally ob- 
served spin-Coulomb, or Pauli spin, blockade^ caused 
by a full occupation of the two-electron triplet state in 
a DQD. Tuning the system such that the level off-set 
Ae — ea — £b = H — £a and the interdot Coulomb re- 
pulsion are half of the on-site charging energy U, pro- 
vide a unit probability for populating the two-electron 
triplet state in one direction of the bias voltage for small 
interdot tunnelling t. The almost full population of the 



two-electron triplet state prevents a current flow through 
the DQD, however, a small leakage current flows through 
the structure because of non- vanishing populations in the 
one-electron states and the lowest two-electron singlet 
state. Depending on the ratio 2t/Ae, the leakage current 
in the blockaded regime is around 1 — 4 orders of magni- 
tude smaller than the current in the near outside of this 
regime. The population of the triplet state is lifted in the 
other bias direction which allows a substantial increase 
of the current through the system. We have also shown 
that a finite ferromagnetic exchange interaction between 
the QDs extends the range of the Pauli spin blockade 
regime, whereas an anti-ferromagnetic exchange yields a 
smaller range of the Pauli spin blockade regime. 
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